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Abstract 

An x-ray pulse-shaping scheme is put forward for imprinting an optical frequency comb onto the fluorescent 
radiation emitted on an atomic x-ray transition, thus producing an x-ray frequency comb. A four-level 
system is used to describe the level structure of N ions driven by small-bandwidth x rays, an optical laser 
field, and an optical frequency comb. By including many-particle constructive interference and the resulting 
enhancement of the emitted resonance fluorescence, a spectrum is predicted consisting of a frequency comb 
centered on the atomic x-ray transition energy and with the same tooth spacing as the driving optical 
frequency comb. This scheme paves the way towards the generation of x-ray frequency combs with energies 
from the 100-eV to the kiloelectronvolt range. 

PACS numbers: 32.50.+d, 32.80.Qk, 42.50.Gy, 42.62.Eh 
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The introduction [1] and development [2] of optical frequency combs led to major advances in 
precision spectroscopy [3-5]. The spectrum of a frequency comb consists of a series of equally spaced 
teeth, namely, the modes of a train of femtosecond pulses spaced by the repetition frequency of a 
mode-locked laser. By counting the number of teeth between a known optical reference frequency 
and an unknown frequency, this comb is used as a very fine ruler to measure an optical frequency 
instead of the corresponding wavelength, which could be determined much less precisely. This 
allows one to reach relative accuracies up to 10~^^ [6]. By precisely counting optical oscillations, 
e.g., in trapped atoms and trapped ions standards, optical frequency combs play a crucial role in 
the realization of all-optical atomic clocks [4, 7]. 

In light of the success of optical frequency-comb metrology, it is highly desirable to extend 
this technology to extreme ultraviolet (XUV) and x-ray frequencies [5]. X-ray frequency combs 
enable precise measurements of high-energy transitions that parallel the accuracy achieved for 
optical frequencies, thus representing an improvement of several orders of magnitude: this allows, 
to name but a few examples, even more stringent experimental tests of quantum electrodynamics 
and astrophysical models [8], and search for variability of the fine-structure constant, to which 
transitions in highly charged ions are predicted to be more sensitive [9]. One can also envision 
ultraprecise x-ray atomic clocks. 

XUV frequency combs have been generated via intracavity high-order harmonic generation 
(HHG) [10]. While in conventional HHG an optical pulse in a gas produces a spectrum containing 
the odd harmonics of the optical frequency, in intracavity HHG a train of coherent optical pulses 
generates an HHG spectrum which in each harmonic line is structured into a fine comb. Based on 
this scheme, Ref. [11] reported the observation of frequency combs at ^ 40 nm 30 eV). The 
required optical peak intensity of ^ 10^^ W/cm^ was obtained with a femtosecond enhancement 
cavity. Yet relativistic effects set a maximum to the highest harmonic for which HHG operates 
efficiently [12]. 

In this Letter we develop a scheme of x-ray pulse shaping, by using an optical frequency comb 
to imprint its structure onto the resonance fluorescence spectrum emitted by an x-ray transition to 
produce an x-ray frequency comb. For this purpose we take advantage of previous investigations 
of many-color schemes of resonance fluorescence in multi-level systems [13, 14]. Examples of x-ray 
pulse shaping are, e.g., studies of electromagnetically induced transparency for x rays [15], for 
which an optical field is used to control x-ray absorption. Anticipated advances in small-bandwidth 
x-ray sources [16-18] will make our approach an attractive alternative, eventually able to reach, at 
relatively low optical-comb peak intensities, energies up to the kiloelectronvolt range. 
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Figure 1. (Color online) (a) An ensemble of ions is driven by x rays (/cx, brown), an auxiliary optical laser 
field (/cl, red), both linearly polarized along the z direction, and an optical frequency comb (/cc, green), 
linearly polarized along the x direction. All fields propagate in the y direction. The resonance fluorescence 
spectrum (/cp, blue) exhibits an induced x-ray frequency-comb structure, (b) Four- level scheme used to 
describe the coherent interaction between He-like ions and the driving fields. 

The geometry used is displayed in Fig. 1(a). An x-ray field 5x(^, ^), an optical auxiliary 
continuous- wave (cw) field S\,{r^ t), and an optical frequency comb SQ,{r^ t), irradiate an ensemble 
of ions. The fields copropagate in the y direction; at time t and position r, they are given by 

e^{v, t) = £q,o{t) COS [Cjqt + (fq{t) + (fq^Q - kq - v] Bq, (l) 

with amplitude £q^o{t)^ carrier frequency cjg, phase (fq{t)^ carrier-envelope phase (pq^o^ wavevector 
kq = {(jjq/c)ey^ linear polarization vector e^, and intensity Iq = \£q^Q\'^ / (Sna) [19]. Furthermore, 
ex = = e^, ec = where e^, and are the unit vectors in the x, ^, and z direction, 
respectively; further, c is the speed of light and a = 1/c is the fine-structure constant. Atomic units 
are used throughout unless otherwise stated. 

The optical bandwidth of 5l(t', t) is so small that it can be entirely neglected. We first assume 
that 5x(^, t) has constant amplitude, £x,o{t) = £x,o^ and constant phase, (px{t) = 0; the effect of 
the x-ray bandwidth is later taken into account through a stochastic approach [20]. Finally, the 
optical frequency comb has constant phase, (fc{t) = 0, and a periodic amplitude [21], £FC,o{t) = 

• 27rfc ^ rji • 27rfc ^ 

J2k=-oG "^k^ 5 with the Fourier coefficients A/^ = 1/Tp ^ ^^fc,o(^) e dt, i.e., 

SFCoii) = ^FCmax En=-oo 0{t - nTp), (2a) 
Git) = C0S2 (t - Si)] R (t _ Si)] , (2b) 

where £^FC,max — ^87ra/FC,max IS its maximum, /FC,max the maximum intensity, and the rectangular 
function is defined with the help of the Heaviside step function 6 as R(x) = 0(x + 1/2) — 0(x — 1/2). 
The full width at half maximum (FWHM) of g^{t) is TpwHM = 2Td arccos ( ^172) /tt [22], with Td 
being the interval in which g(t) is different from 0, <^ Tp. 
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The electric fields 5x(^, t), t), and 5c(^5 are used to drive electric-dipole (El) transi- 

tions in the four-level system depicted in Fig. 1(b), where level i has energy cji and the transition 
energy between the levels i and j is given by Uij — uji— ojj, z, j G {1, 2o,±, 3, 4o,±}. The four-level 
model describes He-like ions, such as Be^^ and Ne^^ [23], with transition energies in the proper 
optical and x-ray ranges. Here the state |1) represents the ground state Is^ ^Sq, with total-angular- 
momentum quantum number J = and positive parity. The three states |2_), |2o), and |2+), with 
Mj equal to —1, 0, and 1, respectively, describe the level ls2p ^Pi, with J — 1 and negative parity. 
Furthermore, |3) represents 1^25 ^Sq, with J = and positive parity, whereas the highest level 
ls2p ^Pi^ with J = 1 and negative parity, is represented by the three states |4_), |4o), and |4+). 
In He-like ions, level 3 has higher energy than level 2 for a nuclear charge Z > 7 [23] [Fig. 1(b)]; 
otherwise, these two levels are inverted in energy. Other levels, such as 1^25 ^^i, ls2p ^Pq, and 
ls2p ^P2, are neglected entirely, since they do not couple via El interaction to the four levels in 
Fig. 1(b) and the spontaneous-decay times from higher-energy levels of Fig. 1(b) to them are by 
orders of magnitude lower than Tp. The excited levels 2, 3, and 4, are nonautoionizing, since in all 
configurations one electron occupies the Is orbital, implying that the levels are energetically below 
the autoionizing threshold [23]. 

The interaction of the three electric fields £q{r, t) [Eq. (1)] with N ions at positions r^, for n G 
{1, . . . , TV}, is described by the Hamiltonian ^ = H^^Y^qH^^q, whereto = Y.n=i I^zg{i,2o,±,3,4o,±} 
is the atomic electron-structure Hamiltonian and H^i^q — 'l2n=i ^n'^qi'^m t) are the El interaction 
Hamiltonians [24, 25]. Here dn = J2i 2o ± 3 4o ±} n-ion, full electric-dipole moment 

operator, with dij = {i\dn\j) being its matrix elements and a^- the n-ion ladder operators. The 
fields are tuned to the respective transition energies, i.e., ux = ^21, = 1^32 1, and uq = ^43. The 
effect of a field on the other transitions to which it is not tuned is entirely negligible: the relevant 
interactions are highlighted in Fig. 1(b). As a result, the undriven states \2±) and |4o), for which 
the decay from higher-energy levels is by orders of magnitude smaller than the decay to the ground 
state |1), can be neglected as well. 

The time evolution of the system obeys the master equation [24] dp/ dt = — i [H, p] where 
p{t) is the density matrix of the system, with elements pfjit) = — Tr{a-^^(t) p}, and the 

Liouville operator C[p] represents the norm-conserving spontaneous decay of the system. 




i,jG{l,2o,3,4±}n=l 

U0i<U0j 



where the El decay rates are given by Fji — Auj'^-^a^\dij\^ /?> [24]; norm-nonconserving terms such 
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as those from autoionization are absent. Further, for cw x-ray and optical fields, and an optical 
frequency comb, the equations of motion (EOMs) are periodic with period T^: the master equation 
admits a periodic solution, p^'^(t) = y5^^(t + Tp), that is asymptotically reached after turn-on effects 
have ceased. 

Though undriven, the El-allowed transition 4^1 undergoes spontaneous decay: because of the 
polarization of the driving fields (Fig. 1), these photons, decaying from states |4±) with Mj = ±1 
to the state |1) with Mj = 0, differ in energy and in polarization from those spontaneously decaying 
on the 2^1 transition. Here we compute the spectrum of resonance fluorescence emitted on the 
4^1 transition with energy CJ41 in the x-ray regime. 

The spectrum of resonance fluorescence is the sum of a coherent and an incoherent part [24, 
26, 27], whose space distribution and intensity are highly affected by many-ion effects. The light 
coherently emitted by many different ions, in fact, adds constructively in the Sy forward direction, 
along which the three driving fields propagate; this implies an enhancement of the spectrum of 
resonance fluorescence proportional to N'^ along the preferential emission direction Sy. Owing to 
the lack of constructive interference, though, the incoherent part of the spectrum does not undergo 
any many-particle coherent enhancement: it is only proportional to N and is distributed in space 
with a dipole radiation pattern displaying very low space directionality [27]. In the forward direction, 
resonance fluorescence is dominated by coherent emission, such that we can neglect the incoherent 
part of the spectrum. The quantum regression theorem and the Wiener-Khintchine theorem relate 
p^^(t) to the spectrum of resonance fluorescence on the 4^1 transition [24, 26, 27], whose coherent 
part is 



2Ke[gll^{t + r) gH^it) e-^-— )n dr dt. 



where gll^{t) = p^^^(t) e-^[^4it-(/cx+/cL+/cc)-rn] ig the relevant element of the slowly varying, n- 
independent, periodic density matrix ^^^(t), and ^41 = d4+i-ecr_ = d4_i-ecr^. Equation (4) contains 
constructive interference of the four allowed paths due to the four combinations Q\^{t + r) gj^i{t)^ 
with j, / G {4_, 4+} [28]. The periodicity of g^^i{t) implies [26] that Sli\^{r Sy , oj) is constituted 
by a frequency comb centered on the x-ray frequency 0041 with the same tooth spacing as the driving 
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optical frequency comb, 

+00 27rm 
5'4i!i4(^ ^y^ ^) = ^ SmS^u- - ^?^) , (5a) 

m=—oo P 

S^ = "-^N^\^j\Z.(t)r-^'il (5b) 

Here, S^-Y^{r Sy ^ uo) duo dA represents the power detected in the energy interval [cj, uo + duo] and in 
a surface element dA = dQ By centered at r = re^. Because of phase matching of the radiation 
scattered by different particles, the photons emitted in the forward direction are focused in a beam 
of mean area AA = 2Ti'^cr^/{uo^iL) [27], where L is the mean length of the ion sample. It follows 
that Pjn — Sm'^Tv'^cr'^ /{uo 41 L) is the power contained in the mth peak of the spectrum. 

We aim at generating an x-ray frequency comb with (a) the same number of peaks, i.e., overall 
width, as the optical comb driving the ions, and (b) with emitted power comparable with that 
of present-day XUV combs generated via HHG [11]. The power in each peak in the comb 
[Eq. (5b)] is proportional to the modulus squared of the mth Fourier coefficient of Q^^iit). The 
properties of a Fourier-series expansion [29] imply that the overall width of the spectrum is inversely 
proportional to the duration of Q^^i{t). This means that, in order to guarantee a wide x-ray comb 
with as many teeth as in the driving optical one, Q^^i{t) must have a pulse shape as short as 
that of the pulses in the optical comb, e.g., by vanishing in their absence. To this extent, we 
introduce the pulse area Q = Jq^^ 1^34+1 ^FC,o(^) dt: in the case of a pulse satisfying the condition 
Q = 2n7r, the system performs an integer number of Rabi cycles [24], such that the population 
and coherences of the highest level in the system are brought back to exactly at the end of the 
pulse [30, 31]. Conversely, for Q ^ 2n7v the interaction with the pulse is followed by a post-pulse- 
exposure spontaneous decay of the highest level. By properly choosing /FC,max to fulfill Q = 2n7r, 
the population of the states |4±) and the off-diagonal terms g^^i(t) Rabi flop only in the presence 
of an optical-comb pulse: the emitted spectrum [Eq. (5)] is given by peaks whose power [Eq. (5b)] is 
the Fourier coefficient of a function which is different from in an interval of duration Tfwhm, and 
has an overall width of the order of 27r/TFWHM- In this way, the long decay Q^^i{t) ^ Q-^^^it ^j^ich 
follows if the intensity is not properly chosen is suppressed: if such an exponential decay takes 
place, in fact, the amplitude of the peaks in the spectrum [Eq. (5b)] is dominated by the Fourier 
coefficients of the function e~^^^^, resulting in a spectrum of smaller width, /41 <C 27r/TFWHM, and 
smaller number of relevant teeth. 

In Fig. 2 we show results from our theory applied to He-like Be^^ ions. The total spontaneous- 
decay rates Fji are calculated with grasp2K [33] and the transition energies UO21 = 121.9 eV, UO23 = 
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Figure 2. (Color online) Evolution in time of the elements of the density matrix g^'^{t) and spectrum of 
resonance fluorescence for Be^+ ions. Present-day parameters are used to model the driving optical frequency 
comb [Eq. (2)], Tfwhm = 120 fs, Tp = 10 ns, 1/Tp = 100 MHz [11, 32], i.e., 27r/Tp = 4.1 x 10"^ eV. The 
ion sample has = 10^ particles over a length L = 1 cm. The driving fields have intensities Ix = 
2.1 X 10^ W/cm^, /l = 1.7 x 10^ W/cm^, and /c,max = 3.0 x 10^° W/cm^, associated with a 27r optical 
frequency-comb pulse. The displayed periodic solutions are (a) g^i{t) for nTp < t < nTp + Td, and (b) 
g^i{t) for nTp < t < {n -\- l)Tp. The power Pm of each peak in the spectrum of Eq. (5) is displayed (c) for 
the whole comb, centered on u^i = 123.7 eV, and (d) around the maximum. In panel (d), ai = 10^ pW~^, 
a2 = 18.64 pW. 



0.2699 eV, 0043 = 2.018 eV, and CJ41 = 123.7 eV, are taken from Ref. [23]. 

By using an optical frequency comb composed of 27v pulses to drive the 3^4 transition, 
population and coherences related to levels |4±) vanish after each pulse. The time evolution of 
^4+1 (^) during a pulse is shown in Fig. 2(a). In the interval in between two optical frequency- 
comb pulses, when the highest states |4±) are completely depopulated and undriven, the remaining 
three states |1), |2o), and |3), behave as a three-level system [13] driven by the fields 5x(^, t) 
and 5l(t', t). During the interaction with a pulse from the optical frequency comb, the elements 
of the density matrix associated with these three states are brought to values, ^^^(T^), which are 
different from the initial ones, ^^^^(0), with z, j G {1, 2o, 3}. To guarantee the overall periodicity 
of the density matrix, ^^^(0) = ^^^(Tp), one takes advantage of the Rabi inversion stimulated by 
the two fields 5x(^, t) and S\,{v^ t), as shown in Fig. 2(b). A proper choice of the intensities /x 
and /l is therefore necessary, in order to maximize the peak of the periodic function ^4^1 (t) and, 
correspondingly, the emitted photon number. 

Having suppressed the post-pulse-exposure spontaneous decay of elements ^4^i(t), the corre- 
sponding spectrum of resonance fluorescence [Eq.(5)] is shown in Fig. 2(c): it is centered on 
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CJ41 = 123.7 eV and contains approximately 10^ peaks of spacing 2ti/T^ = 4.1 x 10~^ eV. Fig- 
ure 2(d) highlights the comb structure of the obtained spectrum. The necessary optical-comb 
peak intensity, /FC,max — 3.0 x 10-*^^ W/cm^, is lower than those presently available and used for 
the generation of XUV frequency combs via HHG [32]. The power of each peak in the emitted 
spectrum — about tens of picowatts — is comparable with that measured in Ref. [11]. 

In order to incorporate the effect of the x-ray bandwidth 7c, we adopt the approach developed 
in Ref. [20]. The x-ray field Sy^{r^ t) is now a stochastic variable, varying in the ensemble of all 
possible realizations: one derives EOMs for the ensemble- averaged density matrix p^^'^^(t) to obtain 
the ensemble-averaged resonance fluorescence spectrum. We find that the 5-peak structure of Eq. (5) 
is now broadened and the spectrum consists of a continuous function whose peaks at ujm — 2T{m/T^ 
have a FWHM of 27c. In order to preserve the comb structure of Eq. (5), the bandwidth of the 
x-ray source must be smaller than the repetition frequency, i.e., 27c < 2ti/T^ — 4.1 x 10~^ eV: the 
many-peak structure would otherwise disappear and the spectrum would reduce to a single wide 
peak. Such a small bandwidth is not available for x-ray frequencies at present. Yet one can weaken 
this condition by increasing the repetition frequency of the optical frequency comb 27r/Tp, though 
this also results in a wider tooth spacing of the x-ray frequency comb. We notice, however, that 
the quality and coherence of x-ray light sources has dramatically improved during the last decades. 
In this respect, the advent of the x-ray free-electron laser (XFEL) [34] has raised encouraging 
prospects, e.g., self- or laser-seeding methods at XFELs [16] or the possibility to use an XFEL to 
drive an inner-shell x-ray laser [17, 18]. Advances in small-bandwidth x-ray sources are, therefore, 
expected to make our scheme a valid alternative for the generation of x-ray frequency combs. 

In conclusion, we present a method for x-ray pulse shaping and use it to generate an x-ray 
frequency comb via the coherent interaction of an optical frequency comb in a multi-level system. 
We have applied our theory to He-like Be^^ at an x-ray energy of 123.7 eV and shown that the 
required peak intensity of the driving optical frequency comb compares with intensities from HHG- 
based methods. The scheme has general validity and can easily find use in modeling different 
physical systems. Qualitatively similar results can be obtained from other He-like ions, such as 
Ne^^; in this case, with uj/^i — 922.0 eV [23], we can predict a comb in the kiloelectronvolt range 
[35], yet for the transition energy 6^43 = 6.679 eV intense optical frequency combs are still not 
available. However, our model can be also applied to different atomic transitions, e.g., Is^ — > Isnp 
with n > 3 in heavier ions, for which experimentally accessible x-ray and optical energies can be 
found; or even to nuclear transitions up to the gamma range. Albeit relying on a driving x-ray 
field of high monochromaticity, our method will take advantage in the near future of the expected 
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advances in better-quality x-ray light sources; x-ray frequency combs will then enable to extend 
all the advantages of optical frequency-comb metrology up to the x-ray regime and allow for the 
development of the next-generation x-ray atomic clock. 

The work of Z.H. was supported by the Alliance Program of the Helmholtz Association 
(HA216/EMMI). 
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